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Abstract. We investigate the pressure effect on the layered Dirac fermion system, 
, which is realized in quasi-two-dimensional organic compound a-(BEDT-TTF)2l3 . The 

' trajectory of the contact points is investigated using the tight-binding model with 

D . the transfer integrals determined by X-ray diffraction experiments. Vanishing of the 

Dirac fermion spectrum, opening of the gap, and pressure dependence of inter-layer 
magnetoresistance are discussed. 
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1. Introduction 

Layered organic conductors, BEDT-TTF [bis(ethyle-nedithiolo)tetrathiofulvalene] 
salts, exhibit various electronic states due to electronic correlation, for example, 
superconductivity, Mott insulator, and charge ordering, under variation of pressure or 
temperature [TJ [2] . 

Recently, Katayama et al. have theoretically suggested that the 3/4-filled a- 
(BEDT-TTF)2l3 salt is a zero-gap system under the uniaxial pressure along the BEDT- 
TTF molecule stack axis (a-axis) [31 H]. In the zero-gap state, the Fermi surface is 
reduced to two points, and the valence and conduction bands contact at two points 
(±fc ), which are not located in highly symmetric points in the two-dimensional (2D) 
Brillouin zone. Those points shift their positions toward the center of the Brillouin zone 
with increasing pressure j5j EJ [7J E] • In the vicinity of each contact point, two bands show 
linear dispersion with tilted cone like shape. In consequence of the cone like dispersion, 
the low-lying excitation properties of the conduction electrons are described by a tilted 
massless Dirac equation [5]. The first principles calculations support this Dirac cone 
structure [QlfTO]. 

The massless Dirac electron system has been also found in graphene which is a single 
layer of graphite [HI [121 [13] . However, in contrast to graphene, a-(BEDT-TTF) 2 i3 is 
not a two dimensional system but a multi-layered bulk material. In a- (BEDT-TTF) 2I3 
crystal, conducting layers of BEDT-TTF molecules and insulating layers of I3 anions 
stack alternatively along the c-axis. Since the conductive layers are separated by the 
insulating layers, interlayer transfer energy is sufficiently small and thus this system has 
a strong 2D nature. 

It has theoretically been explained by Osada [H] that the experimentally observed 
negative interlayer magnetoresistance [15J is due to the zero-mode Landau Level of the 
massless Dirac fermion. In addition, we have proposed that the presence of a tilted and 
anisotropic Dirac cone can be verified using the interlayer magnetoresistance [16] . The 
interlayer magnetoresistance in a- (BEDT-TTF) 2I3 depends on the in-plane magnetic 
field direction because of tilt. 

The hopping parameters of the a- (BEDT-TTF) 2 I 3 system are controlled by 
pressure. Two Dirac cone locations in the Brillouin- zone also move by pressure, therefore 
the distance between valleys would be changed. This feature is very intriguing, because 
the inter- valley scattering effect would be affected by the change of Dirac cone locations. 

In this paper, we investigate the pressure effect on the trajectory of the two Dirac 
points by changing transfer energies within tight-binding model. We also calculate 
the pressure effect on the interlayer magnetoresistance. We use the parameters of the 
Weyl equation estimated by the tight-binding model [3 H] with the transfer integrals 
determined by X-ray diffraction experiments [IT] , and discuss the Dirac fermion merging 
and the gap opening under high pressure region. 

The organization of this paper is as follows. In section 2, we calculate the pressure 
dependence of the Dirac cone parameters in a-(BEDT-TTF) 2 l3 by the tight-binding 
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Figure 1. The conducting layer structure of a-(BEDT-TTF)2i3 , where the unit cell 
is given by the shadowed region. The bonds labeled by al,- • -,a3,bl,- • -,b4 have transfer 
energy. V c (V p ) represents repulsive Coulomb interaction along a(b)-axis 



model. In section 3, we discuss the merging behavior of Dirac points. In section 4, we 
show the exact solution of the Landau level on the tilted Weyl equation. In section 5, 
we calculate the pressure dependence of the interlayer magnetoresistance by using the 
parameters estimated from the tight-binding model. Section 6 gives conclusions of this 
work. 



2. Pressure dependence 

In figure [TJ the basic model describing electronic state in a- (BEDT-TTF) 2I3 is shown 
[3J M, 12 US]. The unit cell consists of four BEDT-TTF molecules named by A, A', 
B and C according to charge disproportionation. To consider the Coulomb interaction 
between molecules, we use the extended Hubbard model which is given by 

+ V ^ a laa a ]/3cT' a j^' a io C a, 
(ia;j /3),<t,(t' 

where % and j denote indices of the unit cell, and a and /3 are indices of BEDT- 
TTF molecules in the unit cell. In the first term, a\ aa (a iaa ) denotes the creation 
(annihilation) operator for the electron of spin a {=[ , \) at the ith site, and t ia .jp is 
the transfer energy between the (i, a) and (j, (3) sites. The second and the last terms 
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denote repulsive Coulomb interactions where U is the on-site interaction, and V a p the 
anisotropic nearest-neighbor interaction. Following Refs. [5J H] , we introduce the effect 
of the uniaxial pressure along the a-axis (P a ) by changing the transfer energies t A 

t A (P a )=t A (0)(l + K A P a ). (2) 

The transfer energies t A and coefficients K A (A = ol, a2, ■ ■ • , 64) are obtained from the 
data at P a = kbar and at Pf, — 2 kbar [17]. We note that the linear variation of 
the transfer energies under uniaxial pressure is reasonable for weak pressure. At high 
pressure, there are deviations from the linear dependence. Since details of pressure 
dependence of the transfer energies, especially at high pressure, are not known, we use 
the linear functional forms for simplicity. Although critical pressures below cannot be 
taken seriously, the purpose of the present study is to demonstrate pressures effects on 
magnetoresistance and merging effects of Dirac fermions. The Coulomb interactions 
U, V a p are treated within the Hartree approximation. The mean field Hamiltonian is 
given by 

^MF = e ^( k ) a k a a a k(3a 



<a,/3> 



(3) 



per 1 S 

where a^aa is the Fourier transform of a ia(7 , (n aa ) is the averaged number of electrons 
(l/N)J2i ( a laa a iaa), & denotes the inversion of the spin a, and d denotes the vector 
connecting the nearest neighbor sites of the unit cell. 
The Hamiltonian ()3]) is diagonalized as 

4 

22e a ^(k)dp r(T (k) = ira{k)d ara {k), (4) 

0=1 

where r is an index of eigenstates with eigenvalues £ rcr , which are arranged by a 
descending order £i CT (fc) > ^2a(k) > ^(fc) > £,4a(k)- d ara (k) is the corresponding 
eigenvector. 

The averaged number of electrons (n aa ) is written by 



i 
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(rW) ? S eM(J(k) )k B T] + 1 ' (5) 
where T is a temperature and ks denotes the Boltzmann constant. The chemical 
potential /x is determined by the condition | = | J2aa ( n av), because of 3/4-filling. 
The parameters U = 0.4 eV, V c = 0.17 eV, and V p = 0.05 eV are chosen [HI [5]. 

In this paper, we consider the zero-gap state. The conduction band (r = 1) and the 
valence band (r = 2) are degenerate at the two points k and — k , and in the vicinity 
of the contact point ko, the Hamiltonian is written by [5] 

H T (k) = r (v%k x + v y k^j a + rv x k x a x + v y k y a y , (6) 
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where k = k — k 0} a x and a y are Pauli matrices and o~q is the identity matrix, r = ± 
denotes the valley index which corresponds with k Q and — k respectively. These contact 
points occur in pairs and can be described by independent degrees of freedom, which 
leads to a twofold valley degeneracy. The valley degree of freedom is called as "valley 
spin". We note that the k x - and k y -ax.es are not along the crystalline a*- and b*-axes, 
respectively (the superscript * means the reciprocal), because the system is rotated in 
order to remove complexity of Hamiltonian. We define the angle made by k x and b* as 

00- 

For convenience, we define the Dirac cone parameters as 



where vq and 4>t represent the tilting magnitude and the direction of the Dirac cone, 
respectively. The parameter a represents the strength of anisotropy coming from non- 
tilting effect. The parameter 7 measures the strength of tilt of the Dirac cone which 
satisfies the relation < 7 < 1 (7 = 1 for non-tilting case). 

Figure [2] shows the pressure dependence of the Dirac cone parameters under the 
uniaxial pressure P a . Figure E]Ja) shows the pressure dependence of the anisotropy 
coming from non-tilting effect. At P a « 8 kbar, the system is almost isotropic (a ~ 1) 
because the hopping parameter t a2 takes almost the same value as Ui and tb2- At 
P a > 10 kbar, a increases with pressure. This growth results in the increase of the 
interlayer magnetoresistance peak with respect to azimuthal angle dependence. In the 
region 4.5 < P a < 8 kbar, a decreases with increasing pressure. Figure EJ^b) shows the 
pressure dependence of the amplitude of Dirac-cone tilting. At P a « 5 kbar, the tilt of 
Dirac cone takes maximum. In the region 5 < P a < 35 kbar, 7 increases with increasing 
pressure, thus, the tilt of Dirac cone decreases. At P a > 35 kbar, the tilt of Dirac cone 
increases again with pressure. Figure [2(c) shows the angle made by k x and crystalline 
b*-axis. In the high pressure region, P a > 20 kbar, this angle becomes almost constant 
and the k x axis is parallel to the crystalline b*-axis. We recall that the pressure P a is 
uniaxial. In the high-pressure region, the transfer integrals t a \ and t a 2 are enhanced 
by the uniaxial pressure, so the energy contour becomes elliptic and shrinks along the 
Figure E^d) shows the azimuthal angle of the tilting direction. In the high 
pressure region, P a > 30 kbar, the tilt of Dirac cone is almost along the fc^-axis. 

3. Merging Dirac points 

As pressure increases, the two contact points approach each other, and then they merge 
into the single point. After merging, the contact points vanish and the gap opens 
between the electron and hole bands. Montambaux et al. have proposed the universal 
2x2 Hamiltonian to describe the motion and merging behavior of Dirac points, and 
they have obtained a semiclassical description of the Landau levels spectrum P,[7],[8]. 
This model describes continuously the Landau level coupling between valleys associated 
with two Dirac points in the vicinity of the merging Dirac points. 
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Here we calculate the trajectory and the gap opening behavior of the Dirac points 
by the 4-band tight-binding model described in section [2j Figure |3] shows the pressure 
dependence of the distance of the Dirac points for several values of U, V c , and V p 
scaled by an interaction parameter x. The contact points exist on wide pressure range. 
However, in some pressure region, they are not located at the Fermi level. Under low 
pressure regions, P a < 3 kbar, the contact points exist but they are not located at the 
Fermi level because of the existence of the hole and electron pockets, which are denoted 
by dashed lines in figure El Under high pressure region, P a > 68 kbar, the contact 
points in the case of the interaction parameter x = 0.8 and x = 1.0, are not located 
at the Fermi level as denoted by dashed line in figure [31 In the case of the interaction 
parameter < x < 0.6, the contact points are located just on the Fermi level in vicinity 
of the critical pressure for the merging of the Dirac points. After merging, the contact 
points vanish and the gap opens between the electron and hole bands. Figure H] shows 
the pressure dependence of the gap between the two subbands. The bandgap depends 
on the pressure linearly. The critical pressure for the merging of the Dirac points 
increases with increasing the parameter x. This critical pressure increase is not general 
behavior for other interaction parameters. In this case, the interactions U, V c and V p are 
taken so that the charge disproportionation pattern becomes the stripe pattern which 
is consistent with the experiment. We do not understand the mechanism of this upward 
shift. But this upward shift suggests that the Dirac fermions are stabilized by increasing 
the interaction parameters. Probably this is associated with the enhancement of charge 
disproportionation. We would like to investigate this point in a future publication. 

4. Exact solution of Landau level on tilted Weyl equation 

As shown in the section [2j Dirac fermions in the a-(BEDT-TTF)2l3 system are described 
by tilted Weyl equation. Reflecting the tilt of the Dirac cone, the Landau level wave 
functions are anisotropic. In this section we derive the exact Landau level wave functions 
of those Dirac fermions under magnetic field. 

First, we rescale the system to remove the anisotropy coming from non-tilting 
effects: v x n x — > vir x , v y 7i y — > vrr y where v = Jv x v y and 7Tj = —ihd Xi + eAi {i = x,y). 
Second, we rotate the system by the angle <pt in the plane so that the tilting direction 
of the Dirac cone would be along with rotated k x axis: 




(8) 



After these transformations, the tilted Weyl Hamiltonian is written as 



H T (p) = vU ] T ((j)t) [rv n' x o- + Tn' x a x + n' y cr y ] U T ((f) t ), 



(9) 




(10) 
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We multiply both sides of the Schroedinger equation Hip = Eip by the operator 
U\((pt) [t^xO'x + 7l 'y a y\ U T ((pt) from the left, and then after some algebra we obtain 

2 



„~.2> 



v 1 — Vq 



E 2 



+ 7T. 



Jl 



IP 



11) 



TV 



We redefine the momentum operator as 



—Ul((p t ) {<J Z + iv a y ) U T 



ip. 



7T, 



7 7T- + - 



vl — vl 



-7T ' 



_/2 



7 



(12) 



-ift 2 / £ 2 , where / 2 



Both ffa and ffy satisfy the commutation relation [tt X} 7r y 
yfh I eB z is the magnetic length. We rewrite equation (fTTj) by the redefined momentum 
operator 

Ik 2 



7 Z 



-TVqC 



l—vt 



+ r 



where E = j-ve. We define the ladder-operator 

L . ... . 



y/2h 



7T, 



ITT 



yj 5 



(13) 



(14) 



which satisfies the commutation relation [a,a^] = 1. In addition, we define the number 



operator N = a)a = 7^2 (7r 2 + 7r^) — §■ We take the eigenstate of the number operator, 
N(p n = ncp n . The eigenstate of Hamiltonian (JSJ) is denoted by n . 

Now we comment on the difference between tilted and non-tilted Dirac cones. If 
the Dirac cone is not tilting, the right hand side of equation ( fl3l) becomes diagonal. 
In the tilted case, the off-diagonal part does not vanish, hence the wave functions are 
linear combinations of <p n and (p n +i- We write the wave function as ip = (u,v) T <p n , and 
substitute this into equation (fT3l) . then we get the relation 

£ 2 - 7 2 r - 7 3 (2n + 1) -T7 2 ?) e- iT9it 

T7 2 £ e ir ^ £ 2 + 7 2 r-7 3 (2n + l) 

For this equation to have a solution, the determinant of the left-hand side matrix must 
be equal zero, so the eigenenergy e becomes 

e = ±v/ 7 3 (2n + 1 ± 1). (16) 

When e = ± a/27%, the wave functions are given by 




0. 



(15) 



Vn 




V2(l + r 7 ) 

When e = ± a/^YTTT) , the wave functions are given by 



(17) 
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V2(l + r 7 ) 
B„ 



-rv e 
t + 7 



+ 




(19) 



a/2(1 + tt) V 

The coefficients A„ and S n satisfy the relation A n e _1< ^ = sgn(n)-B ra which is determined 
from the Schoredinger equation Hif) n = E n ip n . 

Finally, the energy and eigenstate are written as 



E n = sgn (n) \jl hv x v y \ e \ B/f \n\, 



+ 



-TV 

(r + 7 )e iT ^ 



( r + 1 \ 




I -rd e iT ^ ) 


4>\n\-l 



(20) 



(21) 



(n^O), 



and 



r 



-rv e 
^2(1 + r 7 )' I t + ^ 



-lT0t 



(22) 



The explicit form of the Landau level wave functions depends on the choice of the 
gauge. In order to see anisotropy of the wave function, it is convenient to take symmetric 
gauge. On the other hand, for the calculation of the inter-layer magnetoresistance, it is 
convenient to take the Landau gauge. Below we show both cases separately. 



4-1. Symmetric gauge case 

As we shall see later, the interlayer magnetoresistance in a-(BEDT-TTF) 2 l3 depends 
on the in-plane magnetic field direction because of anisotropy in the Landau level wave 
function. In order to get a clear picture, we solve the tilted Weyl equation with the 
symmetric gauge A = x r. 

The presence of the in-plane magnetic field is taken into account by a gauge- 
transformation 



tfi = exp 



x y \ z 



(23) 



with the magnetic length = yfh] eB^ (/i = x,y). After this transformation, the 
vector potential is given by ir x = —\hd x + \eB z y and ix y - 
and rotate the system as 



ihd y — \eB z x. We rescale 



v x 1 x — > vx, v y 1 y — > vy, 



cos 4> t — sin (fi t 
sin (fi t cos <fi t 





(24) 
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respectively. The transformation of equation (1121) is equivalent to x = -j=x' and 
y = y/^ {%)' + -^-^r^f-^- The center coordinate of the cyclotron motion X, Y reads 

X = x- l p y , Y = y+ l ^n x . (25) 

This satisfies the non- vanishing commutation relation as [X, Y] = il 2 , and wave functions 
cannot be simultaneously eigenfunctions of both of them. We choose to use the operator 
X 2 + Y 2 which also commutes with Hamiltonian. In the symmetric gauge case, this 
operator is given by 

-2 , t>2 _ o;2 ' -~ = 1 



X 2 + Y 2 = 2f g \^N-L Z + - (26) 

with the number operator N = a) a , and the angular momentum L z = \{xp y — yp x ) = 
hia^a — tfb). Here, we define the ladder-operator b as 

1 f X + if) . (27) 



y/2l z 

a and b satisfy commutation relations: [b, tf] = 1, [a, b] = [a, b<] = 0. We can rewrite 
these expressions by introducing complex coordinates x — 

a = -^(-ifl e .-iix), l = ^ d x + \r)- (28) 

Using these operators, we find that the eigenstates are denoted by a ket vector 
\n, m) (n > 0, m < —n), where a^a \n, m) = n \n, m) , and tfb \n, m) = (n — m) \n, m). 
The eigenvalue of L z is hm. 

The wave function for zero- mode eigenfunction |0, 0) is obtained by solving a |0, 0) = 
6|0, 0) = 0. In the coordinate representation 4>o,o( r ) = (r* 1 0, 0) , 

MO = ^r e "^ = -4=r ex P (-^jf-) ■ ( 29 ) 



2nl z V2ttI z V 4 ^ 

Higher Landau level wave functions are derived as In, m) = i?LL-_£lp == |0,0). Thus, 

vn! (n—m)\ 

the coordinate representation of wave functions is given by 

<Pn,m(r) = N m , n exp (-^) x lml L^ (^f) , (30) 



4 



where N m n = x / V 1 ; rr is a normalization constant and L™ is Laguerre 

"*>' v27ri z y 2 l - n ~ m > (n— m)\ n 

polynomial I%(t) = ^e*r m |^(e-'r +m ). 

In particular, the n = wave function is written by 

= -=i=-x m exp • (31) 



y/2n2 m m\l z \ 4 

We show |0 o ,o(^)| 2 an d |0o,m(^)| 2 in figure El In the presence of Dirac-cone tilting, the 
energy contour of the cone becomes elliptic. From the uncertainty principle, we expect 
that the wave function shrinks in the tilt direction. In fact, for a tilted Dirac cone, the 
zero energy Landau level wave function is anisotropic and shrinks in the tilt direction 
as shown in figure [51 For a non-tilted Dirac cone, the zero energy Landau level wave 
function is isotropic in real space (not shown). 
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In interlayer magnetoresistance calculation, it is convenient to take the Landau gauge. 



We choose the gauge A a 
and perform gauge transformation 



ByZ + A x ^ , Ay 



-B x z + Ay \ A z 



0, d x A\ 



B : , 



ip = exp 



-12 



y_ 

f 2 

'.r 



11/ . 



Then we obtain 

7T. r 



-ihd x 



zA^\ 



—ihd 7l 



(32) 



(33) 



We rescale and rotate the system as equation (124|) . After this transformation, we take 
—B z y' and A^/ = 0. The transformation of equation (TT21) is equivalent to 



as A^ x ) 



x 



■^=%', V = s/l \ y' + ^^r^JrJ- We choose the operator Y to also commute with 
Hamiltonian. In this gauge, the operator Y is given by Y = l 2 p x /h, the momentum in 
the x direction is conserved in this gauge. Thus, the wave function is given by 

1 ' '~\»ikx ' 



<f>(x,y) 



-4{y)e 



--<P{yy kx 



(34) 



JL X v L x 

where L x is the length of the system in the x direction. The ladder-operator a is given 
by 



V2L 



(35) 



where 



f l = y-Y = ^(y'- l 2 z k) - V2n~v l z . (36) 

The wave function for zero-mode eigenfunction 4>o(y) is obtained by solving a<fro(y) = 0. 
The eigenfunction is given by 



-l) n 



TT^y/2 n n\l z 



exp 



f) 2 
27! 



H„ 



(37) 



where the Hermite polynomials is given by H n (x) 



5. Interlayer Magnetoresistance 

Now we compute the interlayer magnetoresistance and discuss pressure effects on it. We 
represent the magnetic field as B = (B x , B y , B z ) = B(cos 9 cos 0, cos 9 sin 0, sin 9). 

The interlayer tunneling between j z -th plane and (j z + l)-th plane is described by 

H c =-t c £ f d 2 r^ a (ry^ z+1 ^r) + h.c. (38) 

where t c is interlayer transfer energy. In the Landau gauge, the momentum in the y 
direction is conserved because the central coordinate X is conserved, hence the operator 
up is written by 

^( r ) = ^2^n,x(x,y,j z )c x , jz . (39) 

n,X 
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Using equation (1521) . Landau wave function is given by 

1 



exp [izK] ip, ' x (y,z)e 



ipn,x{x,y,z) 



where K represents a phase factor defined by 
K(x, y, 9, 0, ip) = x(sin cos ip — a~ 2 cos sin ip) 



—y[c? sin sin <p + cos cos </?) 
Here the current operator is written by 



J z = iet c — I dx I dy 



h 



B cos e. 



(40) 



(41) 



n,n'k,k' 



e i(k-k')y e ia c K(x,y,9,(j),ip) 



X ^n',k'^n,k€,k,i z d n',k', jz+ 1 + h.C. 



(42) 



where the ^-integration gives k' = k + 5k. Hence, the center of mass Y' is written as 
Y' = Y + l 2 z a^y5k, where 



5fc 



h 



cos 6 1 (sin cos <p — a cos sin (p)a c 



(43) 



The current operator J z becomes 



J z = — / dy exp \ —ia c —B cos 9 (a 2 sin <p simp 

L x I J L ft 

n,n,k,j z '- 

+ cos (j) cos ip)y 

The matrix element (n,k,j z \J z \n',Y',j' z ) is written as 
(n,Y,j z \J z \n,Y',j' z ) 

= iet(5 j+1;jl - o~j-i,j>)o~Y',Y+Pjk 

x J dy exp | —ia c —B cos 9(a 2 sin sin cp + cos cos cp)y 

(X\/y 2 ~ / — 
y - v V2n' ) n ( — j — (y - l z 6k) - v l z V2n 



(44) 



(45) 



t / a V7 
I 



I, 



From the Kubo formula, the interlayer magnetoresistance a zz is given by 

f(E n ,)-f(E n )] \(n,k,j\J z \n'k'j 



EE 

spin i.fc.j 



IUAI\ 1 2 



E n / — -En, J hw + i<5 + (-E n — E n i) ' 



(46) 



where the summations with respect to the layer index j and the wave number k yield 
the number of layer A^ ayer and the Landau level degeneracy = glf, respectively. 
The interlayer magnetoresistance p zz takes the form [16] 

Pzz -Bo 



= 0) 



B + B sin exp 



1 / Oc 

2 I h 



|/(0,a,0 t ,7) 



(47) 
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with 




) 



2 



(48) 



where B is the resistance in the absence of a magnetic field. This formula is derived 
by using the zero-mode Landau level wave function. To justify this approximation, the 
magnetic field B z should be large enough or the temperature is low enough to satisfy 
the relation E\ > ksT. 

The anisotropy of the Landau level wave function shown in figure [5] leads to 
anisotropy in the interlayer magnetoresistance. Figure [6] shows the physical picture of 
the dependence of interlayer magnetoresistance on the in-plane magnetic field direction. 
The in-plane magnetic fields, B x and B y , are treated by the gauge transformation (|32|) . 
which gives rise to the phase factor when the electron hops between one layer to the 
adjacent layer. Figure |6] shows the case that the in-plane magnetic field is parallel to 
the x-axis. In this case, the phase factor is written by exp (iac^j^y). The wave function 
oscillates in real space along the direction perpendicular to the in-plane magnetic field 
because of the phase factor. As a consequence, the matrix element ( l4"5j) is reduced 
when the in-plane magnetic field is perpendicular to the direction in which the wave 
function is extended. Reflecting the real space anisotropy in the wave function, the 
matrix element depends on the direction of the in-plane magnetic field, 0. The inter- 
layer magnetoresistance, thus, depends on 0. 

Figure [7] shows the in-plane magnetic field direction dependence of the interlayer 
magnetoresistance for different pressures. When the pressure increases, the minimum of 
the magnetoresistance moves to 90 degrees and the peak of interlayer magnetoresistance 
increases. In high pressure region, the parameter a increases as shown in figure EJ^a), so 
the effect from the anisotropy coming from non-tilting effect becomes dominant. This 
growth results in the increase of the interlayer magnetoresistance peak. In this case, 
the energy contour shrinks along the a*-axis by the uniaxial pressure, so the interlayer 
magnetoresistance takes the minimum when the in-plane magnetic field is parallel to 
a*-axis, i.e., 90 degrees. 

6. Summary 

In the present study, we examined pressure effects on Dirac fermions in a-(BEDT- 
TTF) 2 l3 within the tight-binding model. The electron and valence bands are degenerate 
at two contact points ko and — fco in the Brillouin zone. They are located at the Fermi 
level under wide pressure range. The pressure dependence of the distance between 
contact points in the Brillouin zone also depends on the interaction parameters. In the 
vicinity of the merging, "valley spin" picture would breakdown because the coupling 
between two valleys, which is usually neglected in graphene, becomes strong. This 
merging behavior may be observed in the pressure range 40 < P a < 70 kbar, and the 
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interaction parameter x < 0.6. Around that pressure, we expect a rapid increase of the 
interlayer resistivity coming from the opening of an energy gap. This suggests that this 
system is useful for investigating valley splitting effect that is still in controversial in 
graphene. We show the exact solution of the Landau level on the tilted Weyl equation 
by using the symmetric and Landau gauges. Because of the tilt, the Landau level 
wave functions become anisotropic and shrink in the tilt direction in real space. We 
calculate the pressure dependence of the interlayer magnetoresistance by using the 
parameter estimated from the tight-binding model. In high pressure region P a > 10 
kbar, anisotropy increases with pressure. This increase results in the increase of the 
interlayer magnetoresistance peak. 
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Figure 2. Pressure dependence of the Dirac cone parameters for U = 0.4 eV, V c = 0.17 
eV, and V p — 0.05 eV. (a) The anisotropy a, (b) the amplitude of tilt 7 (7 = 1 for 
non-tilting case), (c) the angle </>o made by k x and crystalline b* axis, and (d) the 
direction of tilt <pf 
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Figure 3. Pressure dependence of the distance between two contact points. The 
interaction is taken to be U — 0.4a;, V c — 0.17a;, and V p — 0.05x. The dashed line 
denotes that the contact points is not located on the Fermi level. 




Figure 4. Pressure dependence of bandgap. The interaction parameters are taken in 
the same manner as figure [3] 
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Figure 5. Lowest Landau level wave functions |^o,m(^)| 2 
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Figure 6. Physical picture of the dependence of the matrix element on the in-plane 
magnetic field direction. 




Figure 7. The dependence of azimuthal angle dependence magnetic field direction on 
the interlayer magneto resistance for various values of uniaxial pressure. 



